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For a given asymptotically flat initial data set for Einstein equations a new geometric invariant 
is constructed. This invariant measure the departure of the data set from the stationary regime, it 
vanishes if and only if the data is stationary. In vacuum, it can be interpreted as a measure of the 
total amount of radiation contained in the data. 
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Introduction. — The mass of an asymptotically flat 
data for Einstein equations measures the total amount of 
energy contained in the spacetime. The mass is zero if 
and only the spacetime is flat. However, the energy can 
be in different forms: it can be in an stationary regime 
or in a dynamic one. This difference is, of course, phys- 
ically important: gravitational radiation will be present 
only in the second case. The purpose of this article is to 
construct a new quantity that can measure how far the 
data are from the stationary regime. In other words, this 
quantity measures how dynamic are the data and it will 
be zero if and only if the data are stationary. In vacuum, 
the dynamic is produced only by the gravitational field, 
then, in this case, the quantity can be interpreted as a 
measure of the total amount of radiation contained in the 
data. 

The construction is based on a new definition of ap- 
proximate symmetries. An approximate symmetry sat- 
isfies an equation that has always solution for a generic 
spacetime and the solution coincides with the Killing vec- 
tors when the spacetime admits them. The new quantity 
is related to the approximate timelike translation. 

There are many physical situations which are dynamic 
but nevertheless are expected to be described by data 
which are "as close to stationary data as possible" , this 
invariant provides a precise meaning of this idea. An 
important example is the quasicircular orbit for binary 
black holes (see [l| and reference therein). The invari- 
ant constructed here will provide a completely different 
way of finding this class of data and hopefully will help 
to solve the discrepancies between the different current 
approaches. 

Symmetries, approximate symmetries and the con- 
straint map. — Let 5 be a 3-dimensional manifold, and 
let A42,S2,X,C be the spaces of Riemannian metrics, 
symmetric 2-tensors, vectors and scalar functions on S 
respectively. Let hab £ M2, Kab € ^2. The constraint 
map ^ : Ai2 >^ S2 ^ C y. X is defined as follows 

K^-KabK''^\ .,^ 



$ 



hab 
Kab 



R 



-D'K„, 



where Da is the covariant derivative with respect to hab, 
R — h'^^Rab, Rab is the Ricci tensor of hab, K — h°-''Kab 
and a, b, c... denote abstract indices, they are moved with 



the metric hab and its inverse h""^. The set {S, hab, Kab) 
is called a vacuum initial data set for Einstein equations 
ff ^hab, Kab) = on S. 

We can compute the linearization of $ evaluated at 

{hab, Kab) 



D$ 



lab 
Qab 



D'^D'jab - Rab^' - A7 - 
-D'Qab + DaQ - Fa 



H 



and its formal adjoint 



D$* 






DaDb-q - rjRab - ^vhab + Hal 

D^aXb) - D^Xrhab + Fab 



(2) 



(3) 



where 7 = 7afc/l''^ Q = Qabh'^'' and H, Hab, Fa, Fab 
vanished when Kab — (since in this article we will only 
consider the time-symmetric case, the explicit expression 
of these quantities will not be needed). 

The constraint map (^ is important not only because 
it characterizes the initial data for Einstein equations, it 
also gives the Hamiltonian of the theory (see 2] and also 
[3j, (Jl). In particular, the adjoint map ^ gives the right 
hand side of the evolution equations in the Hamiltonian 
formulation (see 5] and reference therein). Moreover, 
£)$* has a remarkably property l6|: the elements of the 
kernel of D^* are the symmetries of the spacetime de- 
termined by the initial data (S, hab, Kab)- That is, if 
(77, X°) satisfies D^*{r],X'^) = then the spacetime wifl 
have a Killing vector <^ and (77, X"") are the projections 
of ^ normal and tangential to the space-like hypersurface 
S, respectively. 

Motivated by this correspondence between the kernel 
of D(f>* and symmetries we introduce the concept of ap- 
proximate symmetries. We will make use of the following 
related operator introduced by Bartnik [3| 



V 



led 
Qpcd 



and its adjoint 
V 



V 

Xa 



= D^ 



1 

Dr 



led 

-DPqpcd 



D$* 



V 

Xa 



(4) 



(5) 



where the dot denotes matrix product. We say that 
(77, X"') is an approximate symmetry if it satisfies the 



equation 



VV*{ti,X°') = 



(6) 



and has the faU off behavior at infinity of the Kilhng 
vectors in flat spacetime. Eq. JB)) is solved for a given 
[hab, Kab)- Note that although the composition D^D^* 
is formally well defined, there is a mismatch in the units 
in the domains and ranges of these operators; that is why 
we use V instead of £)$ in Eq. . 

This definition will be meaningful if: a) Every symme- 
try is also an approximate symmetry b) For generic data 
which admit no symmetry there exist always approxi- 
mate symmetries c) We can uniquely associate each ap- 
proximate symmetry with a symmetry in flat space, that 
is, we have approximate time translation, approximate 
boost, etc. Since every symmetry satisfies D^*{ri, X"-) — 
0, a) trivially follows. In the next section we will prove 
that b)-c) are true for time symmetric initial data: in 
order to simplify the analysis in this article we will as- 
sume that Kab = and X°- ^ at infinity, that is, we 
are only going to consider time translations and boosts 
(in the flat case we have 4 solutions). This is an impor- 
tant special case since it includes all the relevant features 
and difficulties of the equation. The general case will be 
studied in a subsequent article. 

Eq. @ can be derived from a variational principle, it is 
the Eulcr-Lagrange equation of the following functional 



J{rj,X^) 



v*{r,,x-)■v*{v,x-)d^,, 



(7) 



where dfj, is the volume element with respect to hab- For a 
symmetry we have J = 0, for an approximate symmetry 
J > 0. If J, evaluated at an approximate symmetry, is 
finite, we obtain new invariants for the data, which mea- 
sure how far is the approximate symmetry from a sym- 
metry. This will be the case for the approximate boosts. 
However, it turns out, that for the approximate timelike 
translation J diverges. This is because the approximate 
timelike translation grows like r at infinity and not like a 
constant as the timelike translation. The new invariant 
(which will be denoted by A) is precisely the coefficient 
of r in this expansion. 

The variational principle |(7J) is a generalization (which 
includes the lapse function rj and it is of fourth order) of a 
variational principle for X" studied in [g] which is known 
as the minimal distortion gauge (see also the interesting 
discussion in [2|). Finally, we want to point out that the 
operators D^ and Z?$* has been recently used to con- 
struct new kind of solutions for the constraint equations 

[i3II3li3lii 

The time symmetric case. — If Kab = and Xa = o(l) 
then, contracting with {r],X°-) Eq. © and integrating 
by parts, one concludes that D(aXb) ~ 0. Hence X°- is 
a Killing vector of hab- But since there are no Killing 
vectors which goes to zero at infinity (see yy|), we have 



X" = 0. Then, in this case, the operators Q, lO reduce 
to 

Vhab) = DaDb^' - Rab^' - A7, (8) 

r* {V)ab = DaDbV - VRab - ^VKb- (9) 

Using that D''V*{r])ab = -r]DbR/2 and h''''r*{T])ab = 
— 2A?7 — r]R we obtain 

PV*iv) = 2AA77 - D'^D^Rab + vRabR''''+ 

^7^AR + 2RArj+^D''rjDaR (10) 

When R = the last three terms in (|10|) vanish and 
the operator T'T'* has a very simple expression. For an 
arbitrary domain Q C S the boundary term is given by 



l''''V*(T^)abdil^ / V(-iab)vdll+ 

Q. Jn 

{ijDa-/ - jDaV + labD\ - TjD^ab) «" ds, (11) 

an 

where ds is the surface element with respect to hab, and 
ri° is the unit normal of dQ pointing in the outward di- 
rection. 

The fall off conditions on the fields can be conveniently 
written in terms of weighted Sobolev spaces iJL where 
s is a non-negative integer and /? is a real number (see 
ll5j [13 [ill and reference therein, we will use the index 
notation of 0). We say that r] e Hf if 77 e -ff| for aU 
s. The functions in if 2° are smooth in M?" and have the 
fall-off at infinity d^rj = o(r''~ I'l). 

For simplicity, we have not included matter fields in the 
constraint map Q. But then, if the metric is static and 
the topology of the data is E'^ it should be flat. In order to 
have non-flat, vacuum, static metrics (i.e. Schwarzschild) 
we will allow S to have n asymptotic ends, that is, for 
some compact set 51 we have that S \ i~l = '^k=i^k, 
where Sk are open sets diffeomorphic to the complement 
of a closed ball in R'^ . Each set Sk is called an end. We 
will also assume that (5, hab) is asymptotically flat: at 
each end Sk there exists a coordinate system x-l,-^ such 



that we have in these coordinates 



(fc) 



s^, e h:^, 



T < -1/2 



(12) 



where 6ij denotes the flat metric and i,j---, which 
take values 1,2,3, denote coordinates indices. We will 
smoothly extend the coordinate system xjf^, to be zero 
in S\ Sk. 

We say that /3 is exceptional if /3 is a non-positive in- 
teger, and we say that P is nonexceptional if it is not 
exceptional. The manifold {S, hab) will be called static if 
there exist a function v G H^i2 such that V*{v)ab = 0. 
Note that in this definition we allow u to have zeros (hori- 
zons) on S, as, for example, in Schwarzschild initial data. 



We are interested in the kernel of PV* in Ht. Define 
N{p) = dimker {VV* : H^ -^ H°_^) . (13) 

The main result is given by the following theorem. 



Theorem 1. Let (S^hat) be a complete, smooth, asym,p- 
totically flat, Riemannian manifold, with n asym,ptotic 

*/3 



ends. Assume that rj G Hfj satisfies 



VV*T] = 



(14) 



Then, rj £ H'^ . Moreover, assume that (3 is nonexcep- 
tional, then, we have the following: 

(i) IfP<0 then N {(3) =0. 

(ii) // < /? < 1/2 then N{f3) < 1, and N{(3) = 1 if 
and only if {S, hat) is static. In this case we have R ~ 
and the unique static solution v, V*v = 0, has at each 
end Sk the following fall off 



v^^>=o{r-') a'i/ = o(r^-l'') 



(15) 



wht 



ik) 



7^ are a constant and r is given by p2|l . 



(Hi) If I < (3 < 2 then N{P) > An. At each end we 
have the following four linear independent solutions of 
Eq. JUI) 



\k) 



A 



(fc)^(fe) 



^{k), 



vU 



-(fc) 



vU 



(16) 



with ci{k)yfl'(k) ^ -^1/2- Where X(i^\ are constants and 
A(fc) = for some k if and only if {S, hat) is static. 

(iv) In the particular case S* = M'^ (which implies n — 
1) we obtain more information: 

(u')IfO<l3<l then N{f3) < 1 and N{P) = 1 if and 
only if (S, hab) is static and hence flat. 

(nr) Ifl< [3 <2 then N(f3) = 4. 

Proof. The differential operator W* is an elliptic oper- 
ator of fourth order with smooth coefficients. Using the 
decay assumption on the metric p2|l one can easily check 
that PV* is asymptotically homogeneous of degree m with 
m > 4 (for the definition of this concept see for example 
|lfi| . this is the standard assumption on the coefficients 
for elliptic operator on weighted Sobolev spaces, see also 
|l7j). Then, by the weighted Sobolev estimate [Tg], it 
follows that rj e 11% for every s. 

(i)-(ii). We use Eq. (HTJ with 7'^'' = V*{r^Y'' and 
Q, = S to obtain 

I V*{r]T''V*{ri)abdiih= i Ban'^ds, (17) 

Js JdS 

with 



Ba = -2ijDaArj + 2A7jDaV + V*{ri)abD\ ~ -ffDaR, 

(18) 
where the boundary integral is performed in a two sphere 
at infinity at every end. Since d^rj — o(r'^^l'l) we ob- 
tain B = o(r2^-3). If /3 < 1/2 then 2/3 - 3 < -2 and 



the boundary term vanished because ds — O^r"^). We 
conclude that for (3 < 1/2 we have V'P*{ri) = ^=^ 
V*{rj) = 0. That is, for /3 < 1/2 the kernel is not trivial 
if and only if the metric is static. To prove i) we use 
the result that there exists no spacetime Killing vectors 
which go to zero at infinity [13] |lj| • The fall off behavior 
for Killing vectors H15|l was proved in [13. Because the 
constants Vq are non zero, this fall off implies that v 
is unique: assume that there exists another solution v' , 
rescale i/' such that i'' ~ ly — o{r'^ ) at some end, this con- 
tradict H15|l . Finally, we prove that static implies R = 0: 
in 11] it has been proved that static implies that R is 
constant, by our falloff assumption it should be zero. 

(Hi) We use the Fredholm alternative in weighted 
Sobolev spaces (see [l^ and [13 , note that we use differ- 
ent conventions for the weights) to prove the existence of 
the four independent solutions at each end: the equation 



'PV*TJ 

if and only if 



F, with F e H^_ 



will have a solution rj & H 



P 



Fvdn^Q 



(19) 



for all V e H^, such that VV*^ = 0, with (3' ^l~ (3. 

We prove first the existence of r]^ (in the following we 
will suppress the end label {k), all the calculations are 
done in one arbitrary end) . Set 77^ — x^ + 7)^ , then f}^ 
satisfies the equation 



VV*{ff^) = -VV*{x^). 



(20) 



We have that VV*{x^) £ H^_^. Since r < -1/2 we can 
take fj^ G ^1/2 ^^ ^'i- l|2L)(l - From the discussion above 
we have that a solution t)-' G H1/2 ^^ H2U|I will exist if 
and only if the right hand side satisfies the condition 
(|19|l . Since (3 — (3' — 1/2 in this case, we can use ii) 
to conclude that a non trivial v will exist if and only if 
(S, hab) is static. Then, if the metric is not static we 
don't have any restriction and the solutions ij^ exist. If 
the metric is static, we compute 



/ vVV*{x^) = [ x^VV*{v) -f / Ban°-ds, (21) 

Js Js JdS 

= (f Bau'^ds, (22) 

JdS 

where 



Ba = -2iyDaAx^ + 2AvDaX^ + V*{x^)abD''v 

- 2xWaAiy + 2AxWaV, (23) 

and we have used that 'P*{v)ab = 0. Wc use the fall 
off (QS)) for V, the faU off ^^ for the metric and the 
fact that Asx^ = (where A^ is the flat Laplacian) to 
conclude that Ba = o{r~'^), then the boundary integral 
(|22|l vanishes and the solution exists also when the metric 
is static. 



For the solution a we proceed in an analogous way. Let 
A 7^ be an arbitrary constant. We have that VP*{Xr) e 
iJ°_3 (here we use that A^A^r = 0). If the metric is not 
static there exists a solution a £ H^,^ of 



rv*{a) = -vr*{\r). 



(24) 



If the metric is static we can compute condition p9f) as 
we did in in Eqs. ^~ if^ 



I vVr* {\r) = -2 i, j/£>aA(Ar)n" ds (25) 

Js JdS 

(26) 



Ids 



where vq ^ Q is the constant given in H15|l and we have 
used A^r — 2/r. Then, the constant A is zero if and only 
if the metric is static. 

For every end, we have constructed four independent 
solutions, then N{(3) > An. 

(iv) The Fredholni index of W* is given by 



i[P)^N{(3)~N{l-P). 



(27) 



When 5 = M'^ the index t(/3) oiW* is equal to the index 
to(/3) of the flat operator A^A^ (see |l9|). To prove (ii')- 
(iii') we will calculate to(/3) and use t(/3) = /-o(/3)- 

Assume A^A^j/ — 0, that is A^w — Q^ w — /S.su. If 
V & H^ with /3 < 2, then w G -ff^Zj, we use the Liou- 
ville theorem proved in Corollary 1.9 of [l5j to conclude 
that w = i). Then /S.sv = 0, using again this result we 
conclude that for < /3 < 1 we have v — 1, and for 
1 < /? < 2 we have v — 1 , a;-' . Then we have lq {(3) = 
for < /? < 1 and to(/?) = 4 for 1< /3 < 2. 

(ii') We have = io(/3) = i(/3) = N{I3) ~ N{[3'). We 
only need to prove the case 1/2 < /3 < 1. In this case we 
have < /3' < 1/2, using (ii) we get N{(3') < 1, and the 
equality holds if and only if the metric is static. 

(nV) We have 4 = io(/?) = *(/?) = N{(j) - 7V(/3'). and 
/3' < 0, using (i) we get N{l3') = 0. 

D 

The new invariants are given by the constants Xi^^ . 
They can be written as a boundary integral at each end 
Sk 



X(i,\ = - — 



^(fe) 



Svr 



OSk 



n°-DaAa(k) ds 



(28) 



Note that a(fe/) = o{r^^'^) at Sk', for k' ^ k, then if we 
calculate (|28f) for a(fc/) we get zero. Also, if we compute 
this boundary for rj^ we get zero. 

If we assume i? = 0, using Eq. ifTUI) we get another 
representation for A/^,) as a volume integral 



A 



(fc) 






(29) 



Since we have the freedom to rescale a by an arbitrary 
constant, we need to normalize a if we want to compare 



A for different data. In the case of 5* = R^ we have a 
natural normalization. Let hab be a metric on M?. As- 
sume that there exists a smooth family hab{^) such that 
for < e < 1 it satisfies the hypothesis of the previous 
theorem and /lafc(l) = hab, hab{0) = Sab- Then, for every 
e we get a solution a{e). The flat solutions a(0) are the 
constants. We normalize a setting a(0) = 1. With this 
normalization, we can calculate A up to second order in 
e 



A«-^ f R'^'Rabdfi + Oie^) 



(30) 



where R"'' = dRab{e)/deU^o- 

The solutions ry-' provide a coordinate system near in- 
finity. In the non flat case this system is unique up to 
rotations, the translation freedom at infinity is fixed be- 
cause the constants are not solutions of Eq. (|14|l . 

Since V*{ri^) £ -f^-3/9 ^^ have that the functional 
J{r]^) defined in ([7|l is finite. The solutions r]^ are the 
minimum of this functional with the boundary conditions 
ry-' = x^ + o{r^''^) at infinity. The numbers J{ti^) are a 
measure of how far is the metric hab of having a boost 
Killing vector. The fall of conditions of jf are like the 
ones for the boost Killing vectors. In contrast, the so- 
lution a has a different fall-off as the time translation. 
This difference is reflected also in the fact that J{a) is 
infinite, it grows like J{a) ~ Xr at infinity. 
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